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Abstract 

The energy deposition spectra of ^^Ar and ^^Ge in a miniature proportional counter are measured and compared in 
detail to the model response simulated with Geant4. A certain modihcation of the Geant4 code, making it possible to 
trace the deexcitation of atomic shells properly, is suggested. Modified Geant4 is able to reproduce a response of particle 
detectors in detail in the keV energy range. This feature is very important for the laboratory experiments that search 
for massive sterile neutrinos as well as for dark matter searches that employ direct detection of recoil nuclei. This work 
demonstrates the reliability of Geant4 simulation at low energies. 
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1. Introduction 

Nonzero neutrino mass and the existence of dark mat¬ 
ter, whose nature is not known yet, indicate that the Stan¬ 
dard Model of particles is incomplete. A sterile neutrino is 
one of the most natural candidates for dark matter parti¬ 
cles [l| . From this point of view direct laboratory searches 
for sterile neutrinos are of a particular interest. An ex¬ 
periment aiming to search for a hypothetical admixture of 
sterile neutrinos in /3-decay of tritium was proposed Q. 
The admixture can be detected as a specific distortion of 
the tritium electrons energy spectrum in a proportional 
counter. The distortion of interest is expected to be very 
small, therefore it is very important to investigate every 
possible source of a systematic uncertainty. To do this, 
we have simulated the response of a proportional counter 
to the decay of ^^Ar and ^^Ge and compared it to the 
measured spectra. We found that the current version of 
Geant4 Q (10.1 patch 2) used for simulation does not re¬ 
produce all specific features of the energy deposition spec¬ 
tra from those sources. We found also that it was be¬ 
cause Geant4 did not trace accurately a fundamental phe¬ 
nomenon namely the filling of atomic shell vacancies. We 
suggested a code modification (patch) that describes the 
deexcitation of atomic shells correctly. In this work we 
present the measured energy spectra in comparison with 
the model responses before and after applying the patch. 
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2. Measured spectra 

2.1. Decay of^’^Ar and Ge 

To understand the methods of measurement and sim¬ 
ulation, it is necessary to know the radiations that are 
produced after the decay of the two isotopes that were 
employed. ^^Ar and ^^Ge atoms decay by electron cap¬ 
ture followed by deexcitation of the daughters ^^Cl and 
^^Ga. Note that in a gaseous proportional counter only 
the total energy of all Auger electrons is detected with no 
distinguishing of separate Auger-lines. 

^^Ar decays by K-, L-, and M-electron capture. The 
ratios of these modes are L/K = 0.0987 and M/L = 
0.104 0. This leads to absolute fractions K = 0.9017, L 
= 0.0890, and M = 0.0093. Taking into account the X- 
ray branching ratios and the X-ray energies arising from 
vacancies in the various shells (see, for example, 0) leads 
to the left part of the Table [T] If the X-rays and Auger 
electrons deposit all their energy in the counter, three main 
peaks at 2.82, 0.27 and 0.02 keV will appear in a pulse 
height spectrum. Those peaks, named as K-, L- and M- 
peak according to their origin, are used to calibrate an 
energy scale in measurements. 

In case of ^^Ge decay the absolute fractions of K-, L- 
and M-captures are about 0.880, 0.103 and 0.017, respec¬ 
tively Q- Further, taking into account the X-ray branch¬ 
ing ratios and the X-ray energies leads to the right part of 
the Table [TJ Again, if the X-rays and Auger electrons are 
all captured in the counter, three peaks at 10.37, 1.2 and 
0.16 keV will appear there. Note that the “L-peak” is not 
really a single peak but is a structure that contains three 
unresolved lines at 1.299, 1.143, and 1.116 keV. The height 
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Table 1: 

Radiation produced in 

the decay of (left) and 

’^^Ge (right). 

The energies of Auger electrons and X-rays are 

expressed in keV. 

Decay 

Sum energy of 

Energy of 

Percent of 

Sum energy of 

Energy of 

Percent of 

mode 

Auger electrons 

X-ray 

all decays 

Auger electrons 

X-ray 

all decays 

K 

2.823 

0.0 

81.5 

10.367 

0.0 

41.4 

K 

0.202 

2.621 

2.7 

1.143 

9.224 

13.7 

K 

0.201 

2.622 

5.5 

1.116 

9.251 

27.4 

K 

0.007 

2.816 

0.5 

0.107 

10.260 

1.7 

K 




0.103 

10.264 

3.5 

L 

0.270 

0.0 

8.9 

1.299 

0.0 

10.3 

M 

0.018 

0.0 

0.9 

0.160 

0.0 

2.0 


of the two last-mentioned components will vary with the 
probability for capture of the associated X-rays. Thus the 
width of the L-peak, and to some extent its average en¬ 
ergy, depend on the gas composition and the counter di¬ 
mensions. If none of the X-rays are captured, the weighted 
average energy of the L-peak is 1.16 keV. 

2.2. Detector and measurements 

A miniature proportional counter made from a fused 
silica tube, as described in Q and Q, was used for the 
measurements. The tube is a cylinder with length 40 mm, 
inner diameter 4 mm, and wall thickness 1 mm. The inner 
surface of the tube is covered with 1 ^m layer of pyrolithic 
carbon that serves as a cathode. A 12 /rm tungsten anode 
wire is stretched along the cylinder axis. Similar counters 
are used in the SAGE solar neutrino experiment [l^ and 
demonstrate a high sensitivity to X-ray and Auger electron 
emission together with a high long-term stability. 

The energy deposition spectrum after the ^^Ge decay 
has been acquired from the counter filled with a mixture 
of Xe-|-11.0% GeH 4 at the pressure of 703 Torr and the 
avalanche gain of about 10^. The intensity of the ^^Ge por¬ 
tion at the beginning of measurement was about 300 Bq. 
The measurement was done at the SAGE deep under¬ 
ground site in a passive shield The background of 

the counter in this condition was less than 10 counts per 
day above 0.5 keV energy range. The ^^Ar spectrum was 
measured in similar counter filled with pure Xe doped with 
tiny amount of ^^Ar at the pressure of 760 Torr and with 
an avalanche gain of about 5T0^; the intensity was 540 Bq. 
The proportional counter anode was directly connected to 
a charge-sensitive preamplifier with the rise-time of 3 ns 
and decay-time of 30 ^s. The shaping time of the readout 
electronics was about 1 /iS for both measurements. 

Gollected pulse height spectra are shown in Eig.[T] The 
total number of events is about 10® for ^^Ar and about 
1.7T0® for ^^Ge. The threshold in the measurements cor¬ 
responded to 0.1 keV for ^^Ar and 0.2 keV for ^^Ge, so the 
M-peak remains invisible for both isotopes. 


3. Simulation technique 


3.1. Monte Carlo modeling 

The modeling was performed with the Geant4 soft¬ 
ware 0. Geant4 is a Monte Carlo code widely used for 
simulation of particles’ interaction and transport in high 
energy physics, medical physics, space research, and many 
other fields. Geant4 is an open source project, which al¬ 
lows users to easily adopt it for their own needs, e.g. im¬ 
plementing new particles and interactions. If the new fea¬ 
tures are expected to be useful for other users, they can 
be included in the subsequent versions of Geant4. 

The geometry of the real proportional counter described 
above was implemented in a Geant4-based code. The 
counter was filled with a gas corresponding to each mea¬ 
surement, and the ®^Ar or ^^Ge sources were uniformly 
distributed throughout the counter volume. The low en¬ 
ergy interactions were modeled with the Penelope model 
set [Il|, which uses EADL (Evaluated Atomic Data Li¬ 
brary), EEDL (Evaluated Electronics Data Library) and 
EPDL97 (Evaluated Photon Data Library). The follow¬ 
ing processes were modeled: photon interactions (Gomp- 
ton effect, Rayleigh scattering, photo-effect and electron- 
positron pair production), interactions of electrons and 
positrons (bremsstrahlung and ionization) and atomic de¬ 
excitation. The Penelope developers claim the model is 
applicable for energies from 250 eV and higher, for all the 
elements Z = 1 — 100, and the atomic deexcitation works 
only for Z > 5, which is caused by limitations of EADL. 
Therefore, we didn’t take into account the electric field 
because it is practically very hard to simulate low energy 
electron avalanches taking place in the vicinity of the an¬ 
ode. However, the main effect of the electric field presence, 
the so called “end effect”, can be easily taken into account 
by another way, as discussed below. 

Results of Monte Garlo modeling are written to ROOT- 
files [l2|, where the energy deposition inside the counter 
and other values are collected for each decay. One billion 
decays were simulated for each source. 


3.2. Systematic distortion of model response 

A correct comparison of the measured pulse height 
spectrum to the model response is only possible after mod¬ 
eling systematic distortion of the simulated one. The sources 
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Figure 1: The energy deposition spectra of ®^Ar (left) and ^^Ge (right) measured in the proportional counter. 


of systematic distortions for a proportional counter are 
well known. The main systematic factor is a broadening 
of pulse heights caused by the statistical dispersion of the 
number of primary ion-electron pairs. Energy resolution 
R is defined as the full width at half maximum (FWHM) 
of the response divided by the energy E: R = FWHM/E. 
In case of normal distribution its dispersion a and FWHM 
are related as FWHM=2.34cr. The analytic expression of 
the response distortion can be written as 

^OO 

Sj{{E) = / S{x)G{^ = 0,a{x),x — E)dx, 

Jo 

where S{E) is the original normalized electron energy spec¬ 
trum, a, E) is the normal (gaussian) distribution with 
a depending on E, and Sr{E) is the final distorted re¬ 
sponse as a result of convolution. In this particular mea¬ 
surement R was about 25% and 13% for the energies of 
2.8 and 10.4 keV, respectively. 

Another significant source of the response distortion is 
the end effect resulting from degraded events due to the 
decrease of the avalanche gain at the counter ends. It 
can be described in terms of the fraction D of degraded 
events. The fraction D is normally proportional to the 
ends’ volume, where the electric field becomes nonuniform, 
divided by the whole counter volume. In this measurement 
D was about 5%. A special design of a counter (e.g., with 
the length increased up to 200 mm) may depress D down 
to ~1%. Hereinafter the distribution of degraded events, 
caused by the end effect, is assumed to have a rectangular 
shape. In this case, the distortion of the response due to 
degraded events has the form 

So{E) = {l-D)S{E)+D f 

Je X 

where So{E) is the response after the distortion. 

At high rates pileup can also drastically distort the 
response. It can be described in terms of the fraction C 


of indistinguishable coincidences of pulses. The fraction C 
is usually proportional to the counting rate multiplied by 
the pulse width. In the case of these measurements, one 
expects C to be of the order of 0.1% or less. Analytically 
this kind of distortion can be expressed as 

Sc{E) = {l-C)S{E)+C [ S{x)S{E - x)dx, 

Jo 

where Sc{E) is the response after the distortion. 

There are another minor sources of uncertainty that 
can also distort the response. For example, the shape of 
the tube or anode wire may be slightly nonuniform, but it 
is not necessary to account for it. 

4. Results of simulation 

The first simulation immediately revealed nonphysical 
peaks in the model response produced by Geant4. Later 
it became clear that the current version of Geant4 does 
not properly trace the filling of shell vacancies, balancing 
nevertheless the total energy release with additional local 
energy deposition. We suggested a code modification that 
gives a way to solve the problem. The model responses 
both before and after the patch compared to the measured 
spectra are presented below. The physical sense of the 
patch is also described there. 

4.I. Response before patch 

The model response of the counter to ^^Ar and ^^Ge be¬ 
fore the patch was applied to the current version of Geant4 
is shown in jagged thin line (red online) in Fig [5] (upper 
panel). The responses were binned in 50 eV and were 
distorted as described above with resolution R=25% at 
2.8 keV, degraded events D=6%, pileup (7=0.01% — for 
37Ar and i?=13.5% at 10.4 keV, D=5%, (7=0.5% — for 
^^Ge (smooth thin line, blue online). Note, that the pa¬ 
rameters were determined a priori as described above in 
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Figure 2: (Color online) The model response of the counter to decay of ^"^Ar (left) and "^^Ge (right) before (upper panel) and after (bottom 
panel) the patch was applied to the current version of Geant4. The solid black line is the measured spectrum; the jagged thin line (red online) 
is the model response; and the smooth thin line (blue online) is the model response after distortion. In addition, the model response smeared 
by the energy resolution R only is shown separately at the bottom panel (black points). 


subsection o The exact values of the parameters were 
obtained from the response after the patch was applied, 
see subsection 14.31 The counts per energy bin correspond 
to measured spectra. 

In comparison to measured spectra the model responses 
reveal a strong discrepancy. First, in the original model re¬ 
sponse just before systematic distortion the K-peak is split 
onto three main lines. The energy resolution of the counter 
is sufficient to observe the splitting in the experimental 
spectrum of ^^Ge if it would exist. In the spectrum of 
^^Ar those lines would remain unresolved. Second, in the 
model response of ^^Ge the intermediate peak at 2.4 keV, 
that is observed in the measured spectrum, has not been 
reproduced. Thus at this stage of simulation it became 
already clear that the current version of Geantd somehow 
proceeds to fill vacancies improperly. Such a conclusion 
was proved, in particular, by the same splitting that was 
observed in the primary radiations of ^^Ar and ^^Ge gen¬ 


erated by Geant4, just before tracing in the counter (are 
not shown in figures). 

4-2. Geant4 code modification 

A modification was done in the class G4UAtomicDe- 
excitation of the Geant4 that simulates the deexcitation of 
an atomic shells. When a vacancy appeared in some bot¬ 
tom atomic shell after an electron capture or an escape due 
to ionization, it cascades to upper shells. If the energy dif¬ 
ference is emitted with an X-ray only, the current version 
of Geant4 will trace this single vacancy moving properly, 
taking into account all known branching ratios. 

In contrast, when an additional vacancy appears due 
to an Auger electron emission, the Geant4 will not trace 
its further transitions. The energy imbalance inevitably 
occuring is compensated by a so-called “local energy de¬ 
position” being added a posteriori and assumed to be non¬ 
ionizing. 
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We have developed a code modification that makes it 
able to trace the transitions of all possible vacancies. 

4-3. Model response after patch 

The model responses to ^^Ar and ^^Ge after the patch 
was applied to the Geantd are represented by jagged thin 
line (red online) in Fig [5] (bottom panel). The designa¬ 
tions and the systematic distortion are the same as for the 
upper panel. In addition, the model responses smeared 
by the energy resolution R only and just before the end 
effect D and the pileup C distortions are shown separately 
by black points. We did not perform a full-scale optimiza¬ 
tion searching for the minimum of chi-squared. Instead we 
slightly varied the parameters R, D and C near the values 
estimated in the subsection 13 . 21 together with the intensity 
of model response in order to fit it to the measured one 
“by eye”. 

Even without full-scale fit it is clear that after the patch 
was applied Geant4 reproduces the measured spectra much 
better. For instance, the degraded events distribution, 
that has a slope in logarithmic scale in the ^^Ge mea¬ 
sured spectrum, is sloped also in the model response. It 
differs from the distribution of the events degraded at the 
ends of the counter, that is expected to have a rectangular 
shape. The effect is due to the near-wall decays when a 
primary electron releases a part of its initial energy in the 
wall of counter. In particular, in this case the degraded 
events are generated by 8 keV KLL Auger electron, being 
emitted after the K-capture and partially escaping to the 
wall. If it is gone completely, only the energy of deexci¬ 
tation of two L-vacancies of the daughter gallium (about 
2.4 keV in sum) will remain to be detected. This peak is 
clearly reproduced in the model response together with the 
brightly expressed valley between it and the L-peak. The 
Auger electrons at lower energies, appearing near the wall, 
reveal the same effect. Such a behavior, with the slope de¬ 
pending on the electron energy, is in strong accordance to 
the geometry origin of the near-wall degradation. 

In smeared model response to ^^Ge, the 6 keV peak is 
reproduced by the secondary X-ray lines also. This kind of 
radiation is emitted in the filling of vacancies occuring at 
the L-shell of xenon ionized by the primary radiations from 
^^Ge. The peak is formed by the events where secondary 
xenon X-rays escape the counter, taking the L-shell bind¬ 
ing energy of about 4 keV away. The same peak is also 
observed in the measured spectrum, and is commonly re¬ 
ferred to as the escape peak. 

5. Discussion 

One may conclude that after the patch is applied, the 
model response reproduces all visible specific features of 
the measured spectrum both for ^^Ar and for ^^Ge, even 
without optimizing the systematic distortion parameters. 
Not only K- and L-peak with the nearly correct intensity 
ratios are well reproduced, but also specific features of the 


^^Ge spectrum are reflected there, including the near-wall 
degraded events slope and the intermediate peaks. 

The measured spectra of ^^Ar and ’^^Ge in a propor¬ 
tional counter thus provide a convincing proof of reliability 
of the suggested code modification. The proposed patch 
has been incorporated into Geant4 ver. 10.2, making it 
able to simulate Auger cascades properly. From the point 
of view of initial goal of the work one can be sure that the 
major sources of uncertainty are under almost full control. 

In addition, given the widespread use of Geant4, we 
note that obtained results have a broad range of signifi¬ 
cance. The correct tracing of deexcitation makes it possi¬ 
ble to reproduce the response of any radiation detector at 
low energy range. This feature is very important in labo¬ 
ratory experiments searching for massive sterile neutrinos 
as well as for dark matter directly via recoil nuclei (see, 
for example, [l^ and Q)- This work enables physicists to 
have confidence in Geant4 simulations in the keV energy 
range. 
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